PHENOMENOLOGICAL THEORY OF FRICTION IN THE 
QUASISTATIC LIMIT : COLLECTIVE PINNING AND MEMORY 

EFFECTS 

Lyderic Bocquet 1,2 and Henrik Jeldtoft Jensen 2 
1 Laboratoire de Physique, URA CNRS 1325 
E.N.S. Lyon, 46, allee d'ltalie 
69364 Lyon Cedex 07, France 
2 Department of Mathematics, Imperial College 
180 Queen's Gate, London SW7 2BZ 
United Kingdom 

Abstract 

When a elastic body is moved quasistatically back and forth over a 
surface, the friction of the interface is experimentally observed to circulate 
through a hysteretic loop. The asymptotic behaviour of the hysteresis loop 
is approached exponentially. We describe how this behaviour is connected 
to the collective properties of the elastic instabilities suffered by the elastic 
body as it is displaced quasistatically. We express the length scale of the 
exponential in terms of the elasticity of the surface and the properties of 
the rough substrate. The predicted scaling are confirmed numerically. 



Resume 

Lorsqu'un mouvement quasistatique alternatif est impose a un solide 
elastique en contact avec un support rugueux, il a ete mesure experimen- 
talement que la force de friction decrit une boucle d'hysteresis. De plus, 
l'approche de la force de friction vers sa valeur asymptotique est exponen- 
tielle. Nous decrivons comment ce comportement est lie aux proprietes col- 
lectives des instabilites elastiques subies par le solide lorsqu'il est deplace 
quasistatiquement. La longueur de relaxation de l'exponentielle est ex- 
primee en terme de l'elasticite de la surface et des proprietes du support 
rugueux. Les lois d'echelle predites sont confirmees numeriquement. 
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1 Introduction 



The present paper presents a discussion of the hysteretic behaviour of the fric- 
tion force between two solid bodies. We use a one dimensional model to develop 
a phenomenological analysis of friction in the quasi-static limit (zero velocity, 
i.e., below the depinning threshold). In particular, we emphasize the connection 
between macroscopic friction and the microscopic elastic instabilities occuring 
when an elastic body is pulled over a microscopically rough surface. Guided by 
computer simulations of a one dimensional elastic chain, we construct a coherent 
collective pinning picture of friction. In particular, an equation for the relation- 
ship between the friction force and the number of elastic instabilities is derived. 
We use this equation to obtain the hysteretic relation between displacement and 
friction force. Our model is qualitatively in agreement with experiments on elastic 
friction 0. 

1.1 Some open questions in friction 

It is useful to divide the problem concerning the friction between two solid sur- 
faces into a number of sub-problems. At the microscopic quantum level the forces 
have their origin in the electronic potential across the interface between the two 
surfaces. These forces can be measured by various types of scanning probe mi- 
croscopy and can in some cases be calculated by use of a quantum theory for 
the surface atoms of the solids ||. Microscopically the two surfaces may often 
be considered as ideal periodic structures and the interaction across the surfaces 
will be spatially homogeneous and periodic. 

On the other hand, when two macroscopic surfaces are in contact, the forces 
between the surfaces will be due to the roughness of the surfaces on macroscopic 
length scales. The friction force will originate in the interlocking of a number of 
asperities on the two surfaces. If we can neglect plastic deformations (tear and 
wear) of the asperities as they are forced past each other, we may expect that 
it is possible to reduce part of the pinning problem to the interaction between 
macroscopic elastically deformable asperities. We can then neglect the particular 
quantum nature of the microscopic interaction and attempt to establish an effec- 
tive phenomenological theory of the friction between macroscopic elastic surfaces. 
New controlled experiments on artificially constructed elastic surfaces make it ur- 
gent to understand theoretically this model system of friction |]J . Even this vastly 
reduced problem contains some rather subtle open questions. 

Recently most of the theoretical and numerical works on the pinning prob- 
lems have focused on the dynamical properties of friction in order to obtain 
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the friction-velocity relationship above the depinning transition, i.e., when the 
applied force is greater than the critical pinning force ||, |], |5|. In particular 
dynamical Renormalization Group techniques have been successfully applied in 
the low velocity limit by treating the depinning transition as a dynamical critical 
phenomena ||, [5J . 

On the other hand, the litterature on the properties of pinning below the de- 
pinning transition is mainly restricted to the pionneering work of Larkin and 
Ovchinnikov (LO) in the 70's 0. Their analysis first provided a sketch of the 
collective properties of pinning in order to construct the critical pinning force. 
Their predictions were successfully compared to experimental results on the pin- 
ning of vortex lines in superconductivity This success has led some authors 
to consider the behaviour below the critical threshold as essentially solved by 
LO ||. However, in view of new experimental results on quasistatic friction and 
contact line motion ^ |3[], it has become clear that the LO analysis is incomplete 
on the following important points. Firstly the LO argument does not involve 
hysteretic behaviour. Secondly their approach does not take any mechanical in- 
stability into account, while it has been shown that friction cannot arise without 
this crucial ingredient (see || 13] and the discussion below). Moreover, new ex- 



periments have shown the existence of memory effects in the approach to the 
critical pinning force, which are characterized by a length scale [p], ||, [E(|. In the 
experiments of Heslot et al. ||10|| , the introduced length accounts for the "age" 



of the pulled system, while in the experiments of Crassous et al. on friction and 
contact line motion [|, this length characterizes the distance needed by the 
system to reach the stationnary state. Therefore it is clear that although the 
LO analysis captures an essential part of the collective mechanism of pinning, 
the link to the above mentioned experimental observations is still lacking. The 
construction of this link is the main object of the present work. 



1.2 The crucial role of elastic instabilities 



As already pointed out by Tomlinson [11], friction arises in an elastic media due 



to the existence of a multiplicity of metastable states. This multiplicity induces 
mechanical instabilities in the system and leads to dissipation and hysteresis 



The instabilities arise whenever the local force-balance equation describing 
the interlocked surfaces becomes multivalued. As the mutual displacement of the 
two bodies is increased, the asperities are forced against each other. The interface 
force increases linearly as a function of the centre of mass displacement. When 
the local force balance becomes unable to sustain the mutual force between the 



3 



asperities, a local rearrangement of the atoms in the interface will take place. 
If the roughness of the interacting surfaces is sufficiently large, the motion of 
the interface atoms takes place in the form of swift jumps from one metastable 
configuration to another. During such an instability the friction force drops 
precipitously by a certain value, for then again to increase linearly upon further 
increase in the relative displacement of the two bodies. 

Let us recall the microscopic features of the elastic instabilities. This is most 



easily done by use of a single degree of freedom picture |T3|, ||. Consider a particle 
at position x elastically coupled to a position X (one may think of X as the 
center of mass of the elastic lattice). We want to follow the motion of the particle 
as it passes over an asperity modeled by a Gaussian peak in the potential energy. 
The energy of this system is given by 

U = ^k(x - X) 2 + A p exp[-(x/R p ) 2 }. (1) 

The static equilibrium of the system is obtained by solving the equation dll/dx = 
for a prescribed X value, i.e., we have to solve the equation 

9 A t 

K(x-X) = ^feM-(x/R P ) 2 }- (2) 

This equation is solved graphically in Fig. [(]. The important point is that the 
equation has a single valued solution for any value of X as long as the slope of the 
force exerted by the Gaussian peak is everywhere smaller than k. This condition 
is 

^ < «• (3) 

When this condition is not fulfilled (as in Fig. []), the equilibrium position of 
the particle x(X) will become a discontinuous function of X : a mechanical 
instability occurs when the system passes over the defect. Accordingly, energy 
is dissipated (e.g. into the rapid degrees of freedom, like phonons) : instabilities 
induce hysteretic behaviour. However, this single particle mechanism does not 
lead to memory effects. Therefore, in a more general case where an elastic medium 
slides over many defects, one still expects the basic mechanism for friction to be 
linked to the occurence of elastic instabilities; on the other hand, memory effects 
are the indication of the collective character of instabilities. 



1.3 Our aims 



Our aim in this paper is to give a mainly qualitative picture of collective pinning 
in the quasi-static limit in order to understand the memory effects. We stress 
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the fact that as in the experiments of ref. [|TJ, we study the system in a regime 
below the depinning transition : the sliding velocity of the two solid bodies strictly 
vanishes and no dynamical effects are expected in our case. In particular the sys- 
tem will be considered at equilibrium under an applied constraint. The velocity 
dependence of the friction force is therefore not the purpose of the present study 

III- 

We will attack the problem along the following line : 

(i) First we perform numerical simulations of a very simple model of the 
system, in order to check if the latter is able to reproduce the experimental results 
of ref. ]TJ. In particular, the following experimental facts should be recovered: 
existence of an hysteresis loop when the system is pulled back and forth over the 
rough surface; existence of a memory length; exponential approach towards 
the stationnary friction force. 

(ii) Then we construct a "microscopic" scenario of friction by analysing the 
microscopic behaviour of the system. 

(hi) Finally, on the basis of the numerical results, we propose a qualitative 
phenomenological model for the collective properties of friction in the quasi-static 
limit. This model provides the link between the collective equilibrium and out- 
of- equilibrium properties of the system below the depinning threshold. 



2 General features of the numerical results 
2.1 A simple numerical model 

Let us first introduce the model we use for our numerical experiments. For 
simplicity we model the two interacting elastic surfaces by a mobile deformable 
elastic medium in contact with a stationary undeformable rough surface. This is 
clearly a limitation in comparison to two deformable elastic surfaces. However, 
our findings will justify our anticipation that this simplification is inessential. 

Our model consists of a one dimensional string of particles at positions Xi 
where i = 1, N, coupled together through elastic springs all of the same spring 
constant k. The chain, of length L, is assumed to be periodic. We use the 
equilibrium length a of the elastic springs as our unit of length a = 1. The 
particles interact with a set of pinning centres in the form of randomly positioned 
(at positions x\ where i = 1,...,N P ) repulsive asperities of density n p . All the 
asperities are modelled by the same Gaussian potential peak of amplitude A p and 
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range R p . The potential energy of the system can accordingly be written as 



U = U e l + U p i n 

k L ~ l 

= 2 *52( Xi ~~ Xi+i ~ a ) 2 + 

i=i 

L N p 

£EA> exp[-(x,-x^) 2 /^] 

i=l j=\ 



(4) 



We are interested in experiments where dynamical effects can be neglected. 
For this reason we investigate the total force produced by the asperities when the 
elastic chain is moved quasistatically through the asperities. 

The simulations are performed in the following way. We start from the ideal 
lattice configuration of the elastic chain. We then use molecular dynamics an- 
nealing to relax the chain to the substrate potential. This is done by starting the 
particles out with a random distribution of velocities. The Newtonian equation 
of motion derived from the potential in Eq. (Q) is then integrated by the leap-frog 
algorithm |13| |. While integrating the equation of motion we gradually extract 



kinetic energy by scaling the velocities of the particles. In this way we move the 
system towards one of the metastable configurations of low potential energy. 

When this initial relaxed configuration has been prepared we study the friction 
(or "pinning") force induced when the centre of mass (c.o.m.) of the chain is 
gradually forced through the asperities at zero velocity V = 0. The most accurate 
way of numerically simulating the quasistatic displacement of the c.o.m., X, is 
by performing cycles of shift and relaxation. The elastic chain is displaced as a 
rigid body a small amount dx (dx/a = 10~ 5 in our simulations) by replacing all 
the particle positions by Xi i— > Xi + dx. While the c.o.m. is kept fixed at this new 
c.o.m. position X new = X Q id + dx, the chain is next relaxed to the asperities. The 
c.o.m. is kept fixed by always counteracting the force exerted by the substrate 
by an external force f ex t applied homogeneously to all the particles of the chain. 
We have 

1 L dU ■ 

J ext — ^ Z./ Q ^ 

2=1 1 

where f ext used in the present time step is calculated from the positions of the 
previous molecular dynamics time step. This method allows one to follow as ac- 
curately as desired the motion through the background potential. One could of 
course also move the chain by applying an external force larger than the present 
friction force. This method leads to results qualitatively equivalent to the above 
described method. The draw back of simply applying a constant external force is 
that the force from the substrate fluctuates in space. A constant force will there- 
fore sometimes be much larger than the pinning force. This leads to unwanted 
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acceleration effects and make it difficult to remain in the limit of quasistatic 
motion [15 |. 



2.2 Existence of hysteresis 

As in the experiments of ref. HI], we move the elastic chain back and forth over 
the "rough" surface, according to the numerical algorithm discussed above. For 
each c.o.m. position of the chain, X com , we measure the force acting on the 
particles of the elastic chain. This force will be denoted as the friction force, 
Ff. The whole curve Ff vs. X c _ _ m _ is then averaged over many cycles and over 
different random spatial configurations of pinning centres. Fig. [] shows a typical 
numerical result. Starting from a given point, the friction force reaches after a 
finite distance a plateau value, independent of the c.o.m. position. This plateau 
value is the static friction force : it is the maximum value for an external force 
before the solid body (here, the elastic chain) aquires a non vanishing velocity. 
Then, when the system is moved in the other direction, the same plateau value 
with the opposite sign is reached, but following a different curve in the Ff versus 
X c . , m . plane : an hysteresis loop is performed during a cycle. 

As in the experiments, memory effects are clearly observed, since it takes a 
finite length for the system to reach the stationnary pinning force. This defines 
a memory length. 



2.3 Exponential decay of the pinning force 

As shown in Fig. || the approach towards the plateau value for the friction 
force is exponential . For example, when the system is pulled in the X c _ a _ m _ > 
direction, the friction force can be very well fitted after an initial small transient 
distance, by the following relation : 

F f (X c . om .) =F OQ + (F - F„) exp|- Xcom c ~ X ° | (6) 

where Foo is the plateau value for the force, F is the force measured at a given 
point Xq. This relation defines the memory length, (. In our simulation, the 
latter was obtained to be of the order of a few R p , the range of the defects. The 
size of the initial transient regime was obtained to decrease when the density of 
pinning centres n p increases. 

This exponential approach towards the plateau value of the pinning force is in 
perfect agreement with the experimental results of ref. QTJ] . Therefore, the simple 
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numerical model we use should contain the essential of the underlying physics 
leading to the memory effects and exponential decay. We can now analyse in 
greater details the simulation results to give a microscopic picture leading to 
these results. 



3 A first empirical understanding 
3.1 Hooke's law and instabilities 

While in Fig. || we considered an average of the friction force over many different 
initial states, we now focus our attention on a particular realization of the numer- 
ical experiment. A typical non- averaged plot of the friction force as a function 
of the c.o.m. displacement is shown in Fig. || This figure is characterized by 
a saw-tooth behaviour, which splits up into linear increase of the friction force, 
separated by steep decrease in the friction force. 

(i) The linear behaviour corresponds to the reversible linear (elastic) re- 
sponse of the system, when an external constraint is applied in order to impose 
a given c.o.m. displacement. This response is characterized by an elastic suscep- 
tibility, kl, defined by 

dF f = k l ■ dX c , Q _ m , (7) 

relating the measured infinitesimal change in the friction force, dFf, to the c.o.m. 
displacement dX c , _ m ,, according to a simple Hooke's law. The parameter kl was 



first introduced by Labusch [16], in the context of lattice deformations in crystals 



and will be denoted as the "Labusch parameter" in the following. 

(ii) The discontinuities in the force are the indication of a dramatic irre- 
versible transformation occuring in the system. A look at the microscopic tra- 
jectories shows that these jumps in the force are intimately connected to a large 
common displacement of a significant number of particles (typically of order 10- 
50 over 500). The typical displacement Sxo of each particle of the moving block 
was always of the order of one lattice spacing in all our simulations : 5xq ~ a. 
This indicates that, because of the imposed c.o.m. displacement, the local equi- 
librium state of the system becomes unstable and a finite jump occurs towards 
a new local equilibrium state Hl^| . In other words, these jumps occur because 
the system cannot support elastically the applied external constraint anymore, so 
that a new metastable state has to be found in order to release the stored energy. 
Dissipation occurs during these abrupt transitions. In the following, these jumps 
will be denoted as "elastic instabilities". 



S 



The previous scenario is in fact very reminiscent of the single particle case 
[i.e., one spring over one defect), discussed above f|, [17J. In particular, the 



multistability of the metastable postions seems thus to be recovered in the many- 
particle system. However the characterization of the associated instabilities is 
now much more difficult, because of the collective character of the phenomena, 
as will be shown in section 4. 

Let us assume that the elastic chain is pulled in the X com _ > direction. 
We introduce the spatial frequency of instabilities, u, defined as the number of 
instabilities occuring in the system per unit c.o.m. displacement; and AF inst 
the drop in force occuring in a single instability. It is useful to define it to be 
the absolute value of the drop, so that AF inst > 0. The latter is assumed, of 
course, to depend on the c.o.m. position (or equivalently on the friction force) 
at which the instability takes place. According to the previous scenario, a simple 
differential equation for the (averaged) friction force, Ff, can be written : 

= k l -u (AF mst ) (8) 

Usi-c.o.m. 

which accounts for both the elastic response of the system and the finite change 
in the force during instabilities. The notation (...) means a non-equilibrium 
average over many different distributions of pinning centres for a given c.o.m. 
displacement, or equivalently for a given friction force. We omit the brackets 
for the Labusch parameter kl, since the latter is found in the simulations to be 
independent of the friction force Ff. This means that the elastic susceptibility 
of the system remains constant, when the friction force increases, i.e., when new 
metastable states are explored. This rather astonishing fact will be discussed 
briefly in section 5. 

We can see in Fig. |] that no instability occurs until the force has reached the 
value Ff = and a rather long elastic relaxation takes place up to that point 
: then v = for Ff < 0. In fact, this is to be expected since in our numerical 
procedure (involving an alternative forward/backward displacement of the whole 
system like in the real experiments), the initial state is not an equilibrium state 
but the stationnary state of the system when it is pulled in the opposite direction. 
Therefore, until some stress is effectively supported by the system (i.e., Ff > 0), 
the system relax elastically the external constraint and no instability occurs. For 
Ff > 0, the instability process is turned on. The frequency of instabilities is found 
numerically to reach a plateau value in a very short distance (much smaller than 
the lattice spacing, a), followed by a small decrease towards its stationnary value. 
This small Ff dependence will be omitted in the following and the frequency will 
be considered as roughly constant (for Ff > 0). 



9 



3.2 A phenomenological law and the exponential decay 



The averaged change in force during an instability, (AF inst ) does not remain 
constant when the c.o.m. is displaced. This can be seen for example in Fig || 
where the non-averaged friction force is plotted versus the c.o.m. displacement. 
The change in force during an instability is seen to increase when the system is 
pulled quasi-statically, so that (AF inst ) is expected to increase when the friction 
force Ff increases. We conjecture a simple linear relationship between these 
two quantities : 

(AF inst ) = 6F + a F f (9) 
where SFq and a are two phenomenological parameters 

This "phenomenological law" has been checked in the simulations by plotting 
the averaged change in an instability as a function of the friction force measured 
just before the instability takes place. Numerically, this procedure involves a 
simple algorithm which detects the instabilities. The latter relies on the measure 
of the numerical derivative of the friction force, which exhibits a dramatic change 
during an instability. This rough indicator has been checked to work with a very 
good accuracy. We have then averaged the plot over many different realizations 
of the initial conditions to compute for a given friction force, the corresponding 
averaged change in force (AF inst ). A typical result is plotted in Fig. |6|. Except 
in the large force region, the numerical points can be fitted with a good agree- 
ment by a straight line, thus validating the phenomenological relation in eq. (P). 
Physically, the increase of the force release during an instability is understand- 
able. Roughly, when the external constraint increases, the "susceptibility" of 
the system increases accordingly, so that the drop in force during an instability 
becomes larger. This point will however be studied in more details in the next 
section. The increase of {AFi nst ) for large forces in Fig. |6] could be the indication 
of the onset of large fluctuations arising in the close vincinity of the depinning 
threshold, although "critical" avalanches (i.e., involving the whole system) were 
not observed in our simulations. But this problem requires a specific careful nu- 
merical work as done by Pla and Nori ||14|| . This is not the object of the present 
work. 

Combining the phenomenological relation eq. (^) with eq. (§) leads to a 
closed equation for the friction force : 

=kl-v (AF inst ) (10) 
= (k l — u5Fq) — va Ff 

This equation predicts an exponential relaxation of the friction force towards 
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a stationnary value given by 



FT = ^- 6 J± (ii) 

J vet a 

The relaxation length, (, is related to the phenomenological parameter a, through 
the simple relation : 

C = l/i/a (12) 

This exponential behaviour is in agreement with both the numerical results (see 
Fig. |) and the experimental results of Crassous et al. [Q. 

Therefore, the phenomenological relation eq. (|9|) provides an initial "empiri- 
cal" understanding of the underlying physics of friction. However, how this simple 
law emerges from the microscopic picture remains to be clarified. Moreover, the 



results obtained within this approach for the stationnary force (eq. (|TT|)) and 



relaxation length (eq. ([12])) are only useful, if some prediction can be made for 
their dependence on the physical parameters of the system ( e.g. density of pin- 
ning centres n p , strength and range of the pinning centres, A p ,R p ). This is the 
object of the next section. 



4 Microscopic picture : towards collective pin- 
ning 

4.1 Description of instabilities 

In paragraph 3.1 we described the instabilities as an instantaneous collective mo- 
tion of an important number of particles of the system, while the other remained 
more or less stationary. Intuitively, it would be appealing to relate the drop in the 
friction force to the number of jumping particles during the instability. To this 
end, we have studied numerically the microscopic behaviour of the system when 
an instability takes place. Especially, after isolating an instability, we separated 
the system into the block of "jumping particles" and the remaining particles. For 
both parts of the system, we computed the change in the force due to the defects. 
This operation was repeated for a number of instabilities of different amplitudes, 
taking place at different c.o.m. positions. The general "rule" we obtained was the 
following : (i) the change in force measured for the jumping particles only was 
found to be very small; (ii) the global drop in the total friction force observed 
in an instability is caused by the remaining particles, which didn't undergo a 
large displacement. This can be interpreted as follows. When each of the Nj ump 
jumping particles move forward by a large, finite amount d~xo, the c.o.m. of the 
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remaining particles will move backwardby an amount Ax c , , m . ~ — Nj ump 5xq/N, 
in order to keep the c.o.m. of the complete system constant. The corresponding 
drop in force is then simply given by the relaxation of the corresponding elastic 
constraint, according to Hooke's law (see eq. (|7|)) : 

F'after f before = AFi ns i Ax c m 

— _ K L SXQ AT /-| 0\ 

N ly jump v ±tJ / 

Eq. (|1^) provides the desired link between the drop in the force occuring in an 
instability and the corresponding "length" of the instability, i.e., the number of 
particles that jumped, Nj ump . 



4.2 Existence of a coherence length 

The afore-mentionned point (i) - that the change in force measured for the jump- 
ing particles only is extremely small - is very striking. In fact, this can be con- 
sidered as a first indication of the collective character of the instabilities. To 
understand this subtle point, we first recall how an instability occur in the single 
particle discussed by Nozieres and Caroli [L3|]. The system involves 

then only a single particle attached to a spring and interacting with a defect. The 
extremity of the spring is moved adiabatically. For a sufficiently soft spring (see 
section 1.2 for details), the equilibrium position of the particle becomes multi- 
valued over a given range of system positions, i.e., the system becomes unstable. 
The evolution of the particle thus exhibits three phases when it encounters the de- 
fect : first an adiabatic move, corresponding to the elastic response of the defect; 
then a large finite jump towards a point where the reaction force is very small 
: this corresponds to the elastic "instability"; finally, an adiabatic move again, 
where only the tail of the defect's potential is felt by the particle. Accordingly, 
the drop in the force Af during the instability is of the order of the threshold 
force before the jump. Now consider Nj ump particles jumping according to such 
an "individual scenario", the global change in force for the whole system would 
simply be AFi nst ~ Nj ump Af, which is proportional to the number of particles 
involved in the instability. This assumes that the drop in the force originates 
from the particles involved in the instability. But this differs dramatically 
from our observations in the numerical experiments. Indeed according to the 
previous points (i) and (ii), the crucial difference is that in our simulated sys- 
tems, the drop in the force during an instability is only induced by the release of 
the elastic response of the particles not involved in the instability. 

We interprete this point by assuming that the system separates into two 
types of particles : a few particles interact "strongly" with the defects, while all 
the other only interact weakly with the defects, i.e., they can be considered as 
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"strongly correlated" . The equilibrium positions of one (or a few) strongly pinned 
particle may become unstable, the latter thus performing a large finite jump. All 
particles "strongly connected" to this particle then simply follow collectively the 
jump of the first jumping particle. Within this scenario, the change in force due 
to the "jumping part" of the system has no reason to be large. Indeed, most of 
the jumping particles (i.e., except for the one -or the few- strongly "pinned") do 
not interact strongly with the defects, so that their individual change in force 
during their jump is mostly random, justifying the observation facts (i) and (ii). 

There is another indication of the "collective" character of the instabilities. 
Let us come back to the phenomenological relation (Q) for the drop in the force 
during an instability, verified numerically on Fig. |6|. As can be observed in this 
figure, the value at the origin, 5F = AF inst (Ff = 0), is non- vanishing. This 
means that, as soon as instabilities are allowed, a finite (non-zero) number of 
particles are involved in the instability. Therefore, the length scale of the pre- 
viously discussed "strongly connected" particles is non vanishing at equilibrium, 
and is hence an intrinsic property of the system at equilibrium. 

These observations leed us to the notion of "strongly pinned" particles and 
"strongly connected" particles, the latter being characterized by a well defined 
correlation length. In fact, this separation is reminiscent of the analysis intro- 
duced by Larkin and Ovchinnikov (LO) |], [7], of the pinning of vortex lines in 
superconductors. Their analysis is based on the introduction of a "correlation 
length" , characterizing the balance between the elastic energy (of the Abrikosov 
lattice in superconductors) and the pinning energy due to the interaction with 
the defects. This length can be equivalently interpreted as the displacement cor- 
relation length of the system. More precisely, the linear size of the correlated 
volume II is defined as the length scale over which the elastic rigidity of the 
elastic medium is sufficient to counteract the random forces. The relative shear 
induced by the random forces is small within a correlated volume. The particles 
within a correlated volume are strongly correlated in the sense that if a particle 
is displaced a small amount, this displacement will be transmitted by the rigidity 
of the lattice out to distances of the order of the size of the correlated volume. At 
distances larger than the size of the correlated volume the random forces become 
essential. The distortion of the elastic lattice becomes appreciable and particles 
farther apart than are only weakly correlated. In this sense we can think of the 
elastic lattice as being broken into weakly interacting rigid subvolumes. It is at 
the interface between the correlated volumes that the elastic instabilities occur. 

We now recall the (qualitative) LO argument in order to estimate it,. 

In the picture depicted above, the total pinning force on a correlated volume, 
i.e., the force due to the defects, is the sum of many contributions of the same 
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order / ~ A p / R p , but with different signs (since the strongly correlated particles 
only interact weakly with the pinning centres). Therefore, the pinning force 
on a correlated volume involves a statistical summation over all the individual 
contributions and will be of order N^ 2 /o, where N c = n p £l is the number of 
pinning centres included in the correlated volume (we now restrict our study to 
a 1-D problem). Since there are L/£l correlated volumes in the system, the total 
pinning force is found to be of order : 

Fpi n ~ L (^J f (14) 

Since this force only acts for a distance of order R p (the range of the pinning 
potential), before changing randomly, the pinning energy is assumed to behave 
like 



SSpin ~ L (^J fo-Rp (15) 

On the other hand, the increase of elastic energy due to the deformation of the 
lattice inside a correlated volume can be estimated as 5£ e i = \iklnijf) ox 2 , where 
k is the bare spring constant, n = 1/a is the density of particles (a is the spring 
length), and 5x is the distortion distance. The term kjnti, takes into account the 
effective stiffness of ti£l springs contained in a correlated volume. The distortion 
length is expected to be of order of R p , the range of the defect potential. For the 
whole system (involving L/ii correlated volumes), we thus find 

6S a ~L±ka (^) 2 (16) 

Adding eqs. ([15|) and (0), we find the change in energy per unit length associated 
with the randomly distributed defects : 

*H fco (f) 2 +© 1/2/ °^ (17) 

The optimized correlation length is found by assuming that the system evolves 
towards a state which minimizes the costs in energy. This is done by minimizing 
the expression (|T7D with respect to £l, leading to 



2k a R p 

1/2 f 
n P ' to 



2/3 



10 

Replacing (|l"8D into eq. ([14]), we obtain a total pinning force which scales like 

„2/3 f 4/3 

F- ~L- p /0 fl9) 
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Larkin and Ovchinnikov |7| assumed that the pinning force in eq. (|T9| ) is 
the out-of-equilibrium pinning force, defined as the force measured at the 
point where the systems depins, i.e., acquires a non- vanishing velocity. The latter 
corresponds to the symptotic force measured in our numerical experiments (or in 
the real experiments of ref. when the system is pulled quasi-statically. 

But our interpretation of F P i n differs qualitatively from the previous affirma- 
tion. Indeed, the LO argument as depicted here is an equilibrium argument, 
which allows one to find the equilibrium properties of a lattice interacting with 
randomly distributed defects. In this case, the averaged pinning force should van- 
ish (otherwise the system is not at equilibrium). Therefore the "pinning force" 
defined in eq. (fUj) cannot be identified with the out-of-equilibrium pinning force. 
Moreover, the LO argument does not take the elastic instabilities into account, 
while the numerical simulations (and even the real experiments, see [IS] and fl9|| ) 
show that these instabilities play a crucial role in the pinning process. 

In contrast to Larkin and Ovchinnikov, we interpret F P i n of eq. (0), obtained 
along the previous argument, as a typical underlying scale for the fluctuations 
of the reaction force due to the defects at equilibrium : more precisely, it fixes 
the amplitude of the root-mean-square (r.m.s.) fluctuations of the pinning force. 
How this scale is related to the non- vanishing, non-equilibrium pinning force (i.e., 
the plateau value for the force for example in Fig. ||), remains to be understood 
and will be the object of the next sections. 



4.3 Linking the coherence length and the phenomenolog- 
ical relation 

In order to make the connection between the equilibrium properties of the system 
and our empirical understanding of friction depicted in section 3, we make the 
following two conjectures : 

(A) in an instability, the number of particle undergoing a jump is given by 
the length of the correlated volume. In other words, the instability is defined by 
a block-jump of an entire correlated volume. 

(B) when an external constraint is applied, the length of the correlated volume 
increases . 

Let us first discuss these two assumptions. The first one, (A), is intuitively 
understandable. Indeed, according to the picture already discussed in the previ- 
ous paragraph, instabilities are initiated by the few particles which are strongly 
interacting with the defects. When one of these strongly interacting particles 
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becomes unstable, and thus peforms a large jump, all the particles which are 
strongly "attached" to this jumping particle will follow as a whole. 



The second assumption, (B), can be justified in the following way. The corre- 
lation length, £l depends quite strongly on the strength A p of the pinning centres, 
as can be seen in eq. flI8|) (remember that /o ~ Ap/R p ) : the weaker the pinning 
centres are, the longer the correlation length is, because the elastic lattice is then 
less distorted. Now, when an external force is applied, the effective strength of 
the defects will decrease, so that the correlation length is expected to increase 
when the system is moved, justifying conjecture (B). To lowest order in the ex- 
ternal force, we may expect a linear increase of the correlation length, as would 
be obtained by applying linear response theory to the system : 

h(F ext ) = IT + 7 • F ext (20) 

where is the equilibrium correlation length obtained along the LO argument, 
defined in ([18]). F ext is the total external force applied to the system in the 
direction of the motion : F ext = F ext • e, with e the unit vector pointing the 
direction of motion. In our case, the external force applied in order to impose 
a given c.o.m. position simply identifies with the measured friction force, Ff, so 
that Ff = F ext . Note that, at first sight, considerations of symetry would only 
predict a F^ xt dependence of 4 as a function of the external field (because of 
the symetry F ext — > —F ext ). However such a counter- intuitive linear behaviour 
usually occur in degenerate systems, where a careful perturbation theory has 



to be done p0| . This is illustrated for example on the so called "Stark effect" 
for hydrogen, where a first order change in the energy of the first exited state 
(which is fourthfold degenerated) is found when an electric field is applied. In our 
case, the degeneracy stems from the multistability of the metastable states, which 
defines spatial "multipoles" . The linear change in £i may be thought as arising 
from the interaction between these multipoles and the applied field. Moreover, 
the linear guess will be confirmed a posteriori by the consistency of the scenario 
which results from it. 

By use of the two conjectures, we can now understand the "phenomenological 
relation", eq. (U), relating linearly the drop in force measured an instability to 
the friction force. Combining eqs. (|20|) and flilf), and setting Nj ump = n £h(Ff) 
according to the previous argument, we obtain 



AP _ k l n.i L (F f ) Sxp 

inst 



k l n.f£ u 5x ( k l n 7 5x \ 

— n — + — n — r Ff 



(21) 



16 



The proposed scenario thus gives a coherent picture of the mechanism leading 
to the phenomenological law. Moreover, identifying this expression with the 
phenomenological relation eq. ([|), we are left with the following "microscopic" 
expressions for the two phenomenological parameters 5Fq and a : 

c- j-, k l n.fj qu Sxo 

oF = $ 

a = KL y 5x0 (22) 



One of the important consequences of eqs. (p2|) is that it will allow us to 
compute the dependence of the two introduced phenomenological quantities on 
the microscopic parameters of the system (density of pinning centres n p , strength 
and range of the defects A p , R p ). This requires more information on the Labusch 
parameter k^, which will be obtained in the next section. 

For the dimensionless parameter a, the derivation is more subtle since it 
involves the new parameter 7 introduced in eq. (p0[). If we think of this parameter 
in terms of a susceptibility (as usually done in linear response theory), then it is 
determined by the properties of the system at equilibrium [^TJ. Now, according 
to its definition (p0|) , it has the dimensionality of a length divided by a force. But 
for a system of springs interacting with randomly distributed defects, we expect 
the length scale to be fixed by the correlation length scale introduced in the LO 
argument (see eq. flT8|)), while the force should be fixed by a typical value of the 
r.m.s. fluctuations of the pinning force. According to our interpretation of the 
LO argument, this scale is fixed by F pin , introduced in eq. (|TJp. Therefore we 
are left with the following expression for the parameter 7 : 

pequ / pequ \ 1/2 , 

yV ( 23 ) 

Ppin \ n P / L JO 

By replacing eq. ( p3|) into eq. (p2|), one obtains for the phenomenological param- 
eter a the expression 

a ~i nSx °-TT^ < 24 > 

^ J n P 

where the expression for Ei is given in eq. (|TSJ) . This expression can be written 
in a more transparent form : 

a ~ ^ (25) 

"pin 



The predictions of eqs. (|22|) and fl24|) will be compared with the numerical 
results in section 5.3. This will allow us to assess the validity of our approach. 
But we first need to characterize more properly the dependence of the other 
quantities, in particular of the Labusch parameter. This is done now. 
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5 Friction and scaling laws 



5.1 The Labusch parameter 

The Labusch parameter is defined in eq. (|7|) as the elastic susceptibility of the 
system. By definition, this susceptibility measures the response of the pinning 
centres (defects) when the particles interacting through the springs are slightly 
displaced from their equilibrium positions, i.e., when the applied force balances 
the pinning force We performed numerically this "experiment" by displacing 
the c.o.m. by a very small amount and analysed the corresponding change in force 
on each defect. Surprisingly, the change in the reaction force of the defects was 
found to be homogeneously distributed over the defects, thus leading to a picture 
of an individual elastic response of the pinning centres to the diplacement of the 
spring system. 

Therefore, we expect the Labusch parameter to scale like 

kl = n p L ko (26) 

where n p is the density of pinning centres and L the length of the system, ko is 
the individual elastic response of a defect, and can be roughly estimated as the 
probability for a defect to interact with a particle, of order ~ n R p (n being the 
particle density and R p the range of the potential), multiplied by a typical value 
of the second derivative of the pinning potential, of order A p jR 2 v . This leads to 

k ~nR p ^ (27) 

Combining, eqs. (|26| ) and (p7|), we obtain that the Labusch parameter should 
scale like 

k l ~ n p A p R- 1 (28) 

This scaling relation has been checked in the simulations, by varying n p for a 
given set of potential parameters A p , R p , and varying A p , R p for a given density. 
The numerical results are obtained to be in reasonable agreement with these 
predictions. The dependence on n p is shown on Fig. || : the numerical value 
of kl is plotted as triangles; the underlying dotted curve is a straight line with 
slope 1. We explored numerically the dependence on the strength and range of 
the defect potential too. This was done by varying slowly A p and R p for a given 
density. The presented results do not involve however an average over different 
random spatial distributions of the pinning centres. The measured points for 
kl are plotted as circles in Fig. || for the dependence on A p and in Fig. || for 
the dependence on R p : The dotted lines are a guide for the eye to indicate the 
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predicted slopes : slope 1 for the dependence on A p (in Fig. |6]) and slope — 1 
for the dependence on R p (in Fig. |6|). The trend is seen to be correct for both 
dependences. However a more extensive numerical study is still needed for a 
whole set of parameters A p and R p , using different densities of pinning centres. 

It is interesting to note that the picture of an individual response of the 
pinning centres for the elastic susceptibility of the system is coherent with the 
fact that the Labusch parameter is found in the simulations to be independent of 
the friction force (i.e., of the applied constraint). We confirmed this result in all 
our simulations, by checking that the slope of the linear part of the non-average 
friction force, as in Fig. || do not change when the c.o.m. of the system is 
displaced. Apart from the previous qualitative argument, based on an individual 
response of the pinning centres, this result remains up to now quite obscure to 
us. 



5.2 The Frequency of instabilities 

Finally, to complete our phenomenological description, we need to characterize 
more precisely the frequency of instabilities, v. Up to now we do not have a full 
understanding of the mechanism of creation of instabilities, but some predictions 
can however be made. 

In the stationary regime, the following relation is obtained from eqs. (|5p and 

(ID 

k l n£ L (F?°) 5x , s 

k l = v AF mst = v LK N f } (29) 

where FJ° is the plateau value of the friction force. This equality imposes the 
scaling of the asymptotic frequency v to be 

v ~ w (30) 

Note that here we are only interested in the scaling properties of the frequency u, 
so that we forget any friction force dependence of £l to focus on its dependence 
on density of defects and other microscopic parameters : we simply use the fact 
that the stationary value of the correlation length £l(Fj°) has the same scaling 
on n p , etc... as t^ u defined in eq. Ql8|) . We cannot a priori apply the expression 



for v in eq. ( p0|) to the transient regime. 

However an alternative argument can be given, which estimates v from the 
fact that the system is broken up into the weakly and strongly pinned particles. 
Recall the mechanism producing the instabilities in the single particle problem 
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(see 1.2.) : a particle P (position x) is elastically coupled to a position X and 
interact with a defect. In the elastic region, we may approximate the force due to 
the defect by a Hooke's law, with a stiffness kd ~ A p /R p . Then, if the position X 
is displaced by a small amount 8X, the equilibrium position of the particle P will 
be displaced by an distance 5x ~ k 5X/(k + kd), where k is the spring constant. 
If the stiffness of the spring is much weaker than kd (which is a condition for the 
existence of instabilities, see e.g. @), we obtain 5x ~ k bXjkd- Now in our many 
body system, we can identify the position X with the c.o.m. position X c o m , and 
the stiffness of the spring k with the effective stiffness of the correlated volume, 
k ~ kjinti). Moreover the particle P can be identified with one strongly pinned 
particle, since in our microscopic scenario the latter are expected to induce the 
instabilities. An instability will occur when the latter will move over a distance 
of order of the range of the potential : 5x ~ R p . This is equivalent to a change in 
c.o.m. position SX given by 5X ~ n£ L k d R p /k, according to our single particle 
discussion. This argument thus predicts a frequency scaling like 

k 1 

k d nR p f£ u 31 

which is consistent -though not strictly identical- with our first guess, eq. (|30|). 

In particular, this argument predicts a dependence of the frequency of insta- 
bilities on the density of pinning centres as v ~ n^ 3 (see eq. (|18D). The numerical 
results are consistent with this scaling as shown in Fig. |5] : the circles are the 
numerical points; the dotted line has a slope 1/3. 



5.3 Summary of the predictions 

At this stage, we can check the predictions for the phenomenological parameters 
and measured quantities obtained within our scenario : in particular for the two 
phenomenological paramaters 5F and a, for the Labsuch parameter, the plateau 
value of the friction force Fj° and the memory length (. Apart from the pinning 
force and the Labusch parameter, we have mainly focussed our numerical study 
on the dependence of the measured quantities on the density of pinning centres 
n p . Since the predicted scalings on n p are not obvious (see the different powers 
of n p in the equations below !), the comparison with the numerical results as a 
function of n p are therefore expected to be already a drastic test for our scenario. 
But clearly more extensive numerical work needs to be done to explore the whole 
parameter space (n p , A p , R p ). 

First we focus on the phenomenological parameters, 5Fq and a. Their "mi- 
croscopic" expressions are given in eqs. (p2|) and (0). Combining the predicted 
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scaling of the Labusch parameter, as displayed in eq. fl26| ) with the expression of 
the correlation length £ e ^ u , given in eq. (|18D, we obtain the following dependence 
CIS db function of the density of pinning centres, n p : 

a ~ n° p (32) 

In particular, this shows that the slope a of the phenomenological law eq. (|9]) does 
not depend on the density of pinning centres ! Both predictions of eq. (^) have 
been checked numerically. The results are plotted in Fig. |^ : symbols represent 
the numerical points and the dashed lines, the predicted scalings. The agreement 
is seen to be correct for both parameters. This confirmation is crucial, since it 
shows that the proposed scenario for the underlying mechanism of friction (see 
section 4) is coherent. In particular, this justifies our two conjectures (A) and 
(B) of section 4.3. 

Now we turn our attention to the other measured quantities. According to 
eq. (|i~l"D , the stationary friction force is the sum of two term, K^jva and 
5F /a. However combining the "microscopic" expressions obtained for all the 
different parameters involved, eqs. (|22|) , (p4|) , (f25|) , (p0|) , we obtain that the sta- 
tionnary friction force simply scales like F pin , the r.m.s. value of 
the fluctutations in equilibrium. Using eq. (|19D, this leads to 

n 2/3 f 4/3 



nf Af (33) 



On the other hand, using (p4[),(p5[),(|30D, the memory length ( is found to scale 

like 



C =-!- 

VOL 

F 



pin 

~ ~ n 

P V P 



^ Al /3 R- 2 '* (34) 



where the microscopic expressions for F pin and k l have been used (see eqs. 

mm)- 

In fact, an equivalent relaxation length can be introduced, defined as the ratio 
between the asymptotic friction force F^° and the Labusch constant kl : 

poo 

8=^- (35) 
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Since Ff° scales like F P i n , this length exhibits the same scaling as the memory 
length ( : 

§~ C „ n -W A'J 3 Rf* (36) 

On Figs. ^| and |5], we summarize the numerical results for the scalings of 
the different quantities as a function of the density of pinning centres. The 
symbols are the numerical points and the dashed lines indicate the predicted 
scalings. The errorbars were estimated in the following way : for a given (random) 
distribution of pinning centres, the measured quantities were averaged over many 
cycles (forward and backward motion), giving a particular value for the desired 
parameter. Then we averaged again on different random configurations of pinning 
centres : this gives a mean value and a typical errorbar for each quantity. 

For the pinning force and the Labusch parameter, the errorbars are within 
the size of the symbols and thus not displayed explicitly in the figure. For both 
parameters, the numerical results (triangles for and crosses for Fj°) are in 
good agreement with the predicted scalings of eqs. (^) and (|2~8f) : these are 
illustrated by the dotted lines in the figure, with slope 1 for kl and slope 4/3 for 

1700 

b s ■ 

Note however the large errorbars involved in the estimation of the memory 
length, £. Indeed, in contrast to the other parameters, the estimated value of ( 
fluctuates a lot when different configurations of pinning centres are considered. 
Good statistics would therefore involve an average over many different (randomly 
distributed) configurations of defects. Such a procedure requires much more 
numerical work and goes beyond the purpose of this first stage study, devoted 
to construct a coherent picture of friction. The numerical results for £ (circles) 
are however coherent with the proposed scaling (0), as shown in Fig. |6]. The 
underlying dotted line has a slope —1/3, as predicted by eq. ([34]) . 

On the other hand, the estimate of the other length 5 is obtained to be more 
precise, since the errorbars involved in the numerical calculations of Fj° and k,l 
are both quite small. The numerical results (squares) are seen to follow quite 
closely the predicted scaling of eq. (|36D , illustrated in the figure by a line with a 
slope —1/3. 

Finally we checked the dependence of the Labusch paramater and of the 
pinning force as a function of the strength A p and range R p of the defects, as 
given by eqs. ( p6|) and ([33]). The results are summarized in Figs. | and The 
dependence of k l has already been discussed in section 5.1. The results for the 
friction force Fj° are plotted as triangles in both figures. The predictions of eq. 
(|3~3"D are illustrated as dotted lines of slopes 4/3 for the dependence on A p and 
—5/3 for the dependence on R p . Here again, the numerical results are obtained to 
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be in good agreement with the theoretical predictions. Note that the agreement 
of the measured and k l with the predicted scalings on A p and R p implies 
necessarily that 5 follows the predicted dependence on A p and R p as displayed 
on eq. 

Let us emphasize again that we only performed a partial exploration of the 
parameter space {n p , A p , R p }, and a more complete numerical work is still to be 
done. However within the explored window, the trend of all the quantities in- 
volved in the proposed scenario is obtained to follow the theoretical predictions. 
Therefore, the numerical results validate with some confidence the proposed pic- 
ture for the onset of friction. 



6 Discussion and conclusions 



This work focuses on the onset of friction at the interface between an elastic media 
and a "rough" surface. In particular, we study how the friction force evolves 
when the system is pulled quasistatically over the surface. Recent experiments 
have shown that, in this controlled limit, the system exhibits hysteresis and more 
surprinsingly, memory effects These memory effects are characterized by a 
length scale. While it was known for a while that elasticy is able to produce 
hysteresis and dissipation, it is shown here that the memory effects originate 
from the collective character of the induced elastic instabilities. 

The collective character of quasistatic friction is related to the existence of 
a correlation length in the system at equilibrium, which expresses the balance 
between the elasticity of the medium and pinning to the defects. When an exter- 
nal constraint is applied to the system, this length is shown to increase linearly 
with the applied force. This linear increase provides a simple explanation for the 
exponential approach to the stationary static pinning force, as observed numeri- 
cally and experimentally ||]. Moreover the evaluation of this length allowed us to 
make some predictions for the dependence of the stationnary pinning force and 
memory length on the microscopic parameters of the system (density, strength 
and range of the defects of the surface). In particular, it is shown that the value 
of the out-of equilibrium stationary pinning force is fixed by the typical scale for 
the fluctuations of the pinning force in equilibrium. The link between the two 
limits is provided by the analysis of the increase of the correlation length out of 
equilibrium in the spirit of linear response theory. The predicted scalings are in 
agreement with the measured numerical results, validating therefore the proposed 
scenario. 

However, our study only focuses on the quasistatic motion of the elastic media, 
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i. e., when any dynamical effect can be neglected. When the system aquires a finite 
velocity, one may expect that the main lines of the proposed scenario will subsist 
and could provide an interesting alternative route to dynamical friction. Such a 
hope is supported by the analogies between some of our results and the heuristic 
model of Heslot et al. |10[ . They show in particular that experimental results for 



sliding dynamics in the creep controlled regime can be accounted for by a simple 
model, based mainly on the assumption that the properties of the system depend 
on its "age" , 0, for which a "successful" definition was found to be : 



t 

exp 



x(t) - x(t') 

A) 



dt' (37) 



where x is the c.o.m. position of the moving body and D is a "memory length". 
In their work, D is assumed to account for plastic flow, but memory effects 
have been shown in this work to occur even in the elastic limit. This indicates 
that memory effects do subsist above the depinning threshold. However they 
have not been discussed on microscopic grounds up to now. The striking analogy 
between this definition of the age and the exponential decay towards the station- 
nary friction force in our case may be only fortuitous, but this should anyway 
deserve some more careful investigation. However, such a generalization should 
be handled with care. In particular, as already emphasized in paragraph 5.2, our 
understanding of the creation of instabilities is only qualitative. In the dynami- 
cal case (finite sliding velocity), sound waves induced by instablilities may play 
a crucial role. These effects are not included in our discussion of the quasistatic 
limit. 



Finally, it is interesting to note that our approach may be generalized to the 
problem of contact line motion in the quasi-static limit. This problem has been 
investigated experimentely very recently on nanoscales, using the Surface Force 
Apparatus (SFA) technique ||. The experimental results present a striking anal- 
ogy with those obtained in the experiments on dry friction in the quasistatic limit. 
In particular, the force acting on the contact line (related to the contact angle 
of the line) is shown to present hysteresis as a function of the displacement of 
the contact line, defining spatial memory effects characterized by a length scale. 
Moreover, the approach to the stationnary value of the force (i.e., of the advanc- 
ing or receding contact angle) was obtained experimentaly to be exponential too. 
Obviously, all these facts are very similar of those obtained both experimentally 
and numerically for dry friction. In the case of contact line motion, we expect 
the motion to be controlled by the competition between the pinning of the con- 
tact line to the defects of the surface and on the other hand, the "stiffness" of 
the line in the direction perpendicular to the motion. If we neglect the energy 
contribution of the already "wetted" defects, the quasistatic motion of the con- 
tact line reduces to a one dimensional problem and it is thus tempting to apply 
our results. This would lead to a force acting on the contact line, H (denoted 
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as the "hysteresis of the contact line" in the contact line language) which scales 
with the number of defects N p as H ~ N^ 3 (see section 5). This guess is in 
rather good agreement with both numerical results of ref. ||22|| , where a scaling 
H ~ Np 7 was obtained, and even with the experimental results of Di Meglio [ 19] , 
where a scaling H ~ N^ 8 was measured. Moreover in the latter reference, it was 
explicitly stressed that avalanches were playing a crucial role in constructing the 
hysteresis properties. These avalanches are very reminiscent of the instabilities 
observed in our numerical simulations. Beyond this encouraging agreement, fur- 
ther work is needed to understand more deeply these analogies and characterize 
more precisely the regime in which these analogies are relevant. 

More generally, it is interesting to note that the nature of instabilities (here 
elastic instabilities) does not play a crucial role in the qualitative understanding 
of the onset of friction. The hope would then be to be able to include plastic 
effects into the scenario too, along the same lines as elastic effects. Work along 
these lines is in progress. 
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Figure 1 : Graphic resolution of eq. (fj) (in the multistable case 2A p jR 2 v > k). 
The full line is the derivative of the defect gaussian potential. The slope of the 
straight line is k, the stiffness of the spring. The intersection gives the solution of 
eq. (@). Black dots are the stable solutions, the open dot is the unstable solution. 

Figure 2 : Hysteresis of the averaged force as a function of the center of mass 
displacement X com when the elastic chain is pulled back and forth over the defect 
surface. The parameters are : length of the elastic chain N = 500; density of 
pinning centres n p = 0.5; parameters of the gaussian defect potential A p = 0.06, 
R p = 0.25. 

Figure 3 : Log-linear plot of the averaged force as a function of the center of 
mass displacement X com for two different densities (top curve n p = 0.35, bottom 
curve n p = 0.2). The dotted lines are an exponential fit of the averaged force 
in the large X com limit. The other parameters are iV = 500, A v = 0.06 and 
R p = 0.25. 

Figure 4 : Plot of the non-averaged force as a function of the center of mass 
displacement, i.e., for one realization of the cycle. The parameters are N = 500, 
n p = 0.2, A p = 0.06 and R p = 0.25. 

Figure 5 : Average of the discontinuity in force during an instability, (AF inst ) 
as a function of the friction force Ff. The dotted line is linear fit of the numer- 
ical points (the few last ones excepted). The slope gives the phenomenological 
parameter a, while its value at the origin gives 5F Q . For the numerical parame- 
ters under consideration (N = 500, n p = 0.8, A p = 0.06, R p = 0.25), this gives 
a = 0.095, 5F = 1.3. 

Figure 6 : Log-Log plot of the two phenomenological parameters a (squares) 
and SFq (circles), as a function of the density of pinning centres n p . The straight 
lines show the predicted scalings : n° for a and n^ 3 for SF . The length of the 
chain is N = 500 and the parameters for the potential of pinning centres are 
A p = 0.06, R p = 0.25. 

Figure 7 : Log-Log plot of the frequency of instabilities as a function of the 
density of pinning centres n p . The straight line has a slope 1/3. The length of 
the chain is iV = 500 and the parameters for the potential of pinning centres are 
A p = 0.06, R p = 0.25. 

Figure 8 : Log-Log plot of the different measured quantities as a function of 
the density of pinning centres n p . The triangles stand for the Labsuch parameter 
Kl] the crosses for the stationnary static pinning force _Fj?°; the circles for the 
memory length £; and the squares for the length 5. The straight lines indicate 
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the scalings predicted within the proposed scenario. From bottom to top, the 
slopes of the lines are —1/3, —1/3, 2/3, 1. The length of the chain is N = 500 
and the parameters of the potential of pinning centres are A p = 0.06, R p = 0.25. 

Figure 9 : Log-Log plot of the stationnary pinning force (triangles) and of the 
Labusch parameter (circles) as a function of the strength of the pinning centres, 
A p . The straight lines indicate the predicted scalings : Fj° ~ A^ 3 ; k l ~ A l p . 
The length of the chain is N = 500 and the density of pinning centres n p = 0.5. 

Figure 10 : Same as Fig. |6| but for the dependence on the range of the 
potential, R p . The straight lines indicate the predicted scalings : ~ R p 5 ^ 3 ] 
kl ~ Rp 1 - The length of the chain is iV = 500 and the density of pinning centres 
n p = 0.5. 
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